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Using first-principles density functional tlieory we predict a ferroelectric ground state - driven 
by off-centering of tlie magnetic Mn^"*" ion - in perovskite-structure BaMnOs . Our finding is sur- 
prising, since tlie competition between energy-lowering covalent bond formation, and energy- raising 
Coulombic repulsions usually only favors off-centering on the perovskite B-site for non-magnetic 
dP ions. We explain this tendency for ferroelectric off-centering by analyzing the changes in elec- 
tronic structure between the centrosymmetric and polar states, and by calculating the Born effective 
charges; we find anomalously large values for Mn and O consistent with our calculated polarization 
of 12.8 ^C/cm^ Finally, we suggest possible routes by which the perovskite phase may be stabilized 
over the usual hexagonal phase, to enable a practical realization of a single-phase multiferroic. 

PACS numbers: 71.15.Mb, 71.20.-b, 75.47.Lx, 77.80.-e 
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Almost sixty years ago, Matthias observed that fer- 
roelectricity occurs in the ABO3 perovskite structure 
when the i?-site is a transition metal cation with a non- 
magnetic dP electronic structure.^ Since then, many such 
perovskite ferroelectric compounds have been iden- 
tified, and the requirement for "d^-ness" has been ex- 
plained in terms of covalent bond formation between 
empty transition metal d and filled O 2p orbitalsi^i^ This 
Matthias rule is problematic however for the design of 
new multiferroic perovskites, since magnetism is most 
readily accommodated in the perovskite structure by par- 
tial occupation of the transition metal d orbitals. Indeed, 
the contra-indication between i?-site ferroelectricity and 
_B-site magnetism^ has prompted the search for and iden- 
tification of new mechanisms for ferroelectricity, which 
do not involve i?-site cation off-centering: these include 
lone pair stereochemical activity,— spin spirals,— charge 
orderingfi and geometric ferroelectricity In this work 
we discuss instead the circumstances under which the 
Matthias rule can be circumvented and magnetic i3-site 
ions can be made to off-center, using perovskite-structure 
barium manganite as our example. 

We begin this report with a review of the physics of 
the second order Jahn- Teller effect, which explains why a 
non-centrosymmetric distortion - required for ferroelec- 
tricity - is usually only favorable for rf" transition metal 
cations. Then we discuss the balance between ferroelec- 
tric and competing non-ferroelectric distortions, which 
in perovskite oxides are often rotations of the oxygen 
octahedra. The main part of the paper is a detailed 
first-principles study of BaMnOs in the metastable per- 
ovskite structure: We find a ferroelectric ground state, 
with a substantial ferroelectric polarization, driven by 
off-centering of the magnetic Mn'^^ ion. Finally, we dis- 
cuss possibilities for realizing this structure experimen- 
tally 

The tendency of a material to ferroelectric instabil- 
ity can be understood within the framework of vibronic 
coupling theory,^ where it appears as the second-order 
terms in the perturbative expansion of the total energy 
with respect to distortions from a high symmetry ref- 



erence phase. As a result it is often called the second- 
order Jahn- Teller (SOJT) effect. Expanding the Hamil- 
tonian as a function of normal coordinate Q about the 
electronic Hamiltonian for the high symmetry reference 
phase, n'-°\ gives 
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Ti^^^ and TiP'^^ capture the vibronic coupling between the 
displacements of the ions from their positions in the high 
symmetry phase and the electrons. Using standard per- 
turbation theory, the energy can be expanded as a func- 
tion of the normal coordinate about the high symmetry 
reference structureiSiiii with energy as 
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Here |0) is the the lowest energy solution of Ti'"' and 
the |n)s are excited states with energies The first- 

order term, (0| H^^^ |0) Q, describes the regular first-order 
Jahn- Teller theorem. This term is non-zero only for or- 
bitally degenerate states, and in the case of d orbitals 
it always leads to centrosymmetric distortions, therefore 
it does not give rise to ferroelectricity. Note, however, 
that in cases where it is non-zero it dominates over any 
non-centrosymmetric second-order distortions. 

In non-orbitally degenerate systems, competition be- 
tween the two second-order terms, which are of op- 
posite sign, determines whether a non-centrosymmetric 
off-centering is favored or not. The first of the two 
second-order terms describes the short-range repulsive 
forces which would result if the ions were displaced with 
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the electrons frozen in their high symmetry configura- 
tion. Since (0| Ti*^^^ |0) is always positive, it always raises 
the energy of the system, and so polar distortions are 
more likely to be favored if this term is small; this 
tends to be the case for closed-shell dP cations with- 
out valence electrons. ^The second of the second-order 
terms, — ^g(^)-i;"oi Q^i describes the relaxation of 
the electronic system in response to the ionic displace- 
ments through covalent bond formation. It is always 
negative unless it is zero by symmetry, and so it fa- 
vors ferroelectricity when its magnitude is large for non- 
centrosymmetric distortions. This occurs when the prod- 
ucts of the ground and lowest excited states are of odd 
parity, so that the matrix elements are non-zero in the 
cases when i?*^") —E^^^^ is small. For Mott insulators with 
partially filled d shells, the top of the valence band and 
bottom of the conduction band are both composed pri- 
marily of transition metal d states. Therefore the ground 
and low-lying excited states have the same symmetry, 
their product with Ti.'^^^ is odd, and the matrix element 
n) is zero. Conversely, for perovskites, the top 
of the valence band is made up largely of O 2p states, and 
the bottom of the conduction band of transition metal 3d 
states, thus the product of the ground and low-lying ex- 
cited states with H'^^^ is even, and the matrix element 
|n) is non-zero. 

Consequently, the balance between the positive and 
negative second-order terms usually results in off- 
centering for c?" cations, such as Ti*^ in the prototyp- 
ical ferroelectric BaTiOs. Here a strong increase in O 2p 
- Ti 3d hybridization accompanies the distortion from 
the high symmetry to the polar structure^^ and so the 
relevant {0\'H''^^ \n) matrix elements are large. This re- 
arrangement of the electrons through covalent bond for- 
mation leads to Born effective charges, Z* = SP/6u, 
which are significantly larger in magnitude than the for- 
mal charges on the ions; these are referred to as anoma- 
lous Born effective charges (BECs). An anomalous BEC 
is therefore a good indicator of the tendency of an ion 
to off-center and is often taken as a signature of fer- 
roelectricity. The calculated BEC of Ti in the high- 
symmetry cubic phase of BaTiOa, for example, is almost 
+7 whereas the formal charge is This is because, 

as the Ti ion moves towards an oxygen carrying its pos- 
itive charge, negative electronic charge flows towards it 
from the oxygen ion, leading to a larger polarization than 
would arise from the ionic component alone. Conversely 
the balance between the second-order terms usually disfa- 
vors off-centering in transition metals with partially filled 
d shells, since in this case the repulsive Coulomb inter- 
actions are stronger than any energy gain from chemical 
bond formation. We note, however, that second-order 
Jahn- Teller behavior is referred to as an effect, not a the- 
orem, because there is no fundamental requirement that 
the guidelines described above always hold. (In contrast, 
first-order Jahn- Teller is a theorem since it is not a com- 
petition between two terms of opposite sign.) 

In addition to this competition between bond forma- 



tion and Coulomb repulsion that determines the ten- 
dency for off-centering distortions, other centrosymmet- 
ric instabilities compete with ferroelectric distortions in 
determining the ground state. We mentioned already 
first-order Jahn- Teller distortions which always domi- 
nate if they are allowed by symmetry. In addition, 
in perovskite-structure oxides centrosymmetric antifer- 
rodistortive rotations of the oxygen octahedra are com- 
mon. When these rotational instabilities are strong - 
usually in cases when the A-site ionic radius is smaller 
than ideal for the perovskite structure - the tendency 
to form a polar ground state is drastically reduced even 
in materials with strongly anomalous B-site BECs. In 
fact, frustration of these rotational instabilities has been 
proposed as a viable route to novel ferroelectricity and 
multiferroic behavior J^ii^ 

In this manuscript, we use first-principles calculations 
to identify a material - perovskite structure BaMnOa - 
in which a polar mode with transition-metal off-centering 
is the dominant instability in spite of the partially filled 
d"^ manifold on the _B-sitc. Our motivation for choos- 
ing perovskite BaMnOa is two-fold: First, CaMnOa, 
which is paraelectric down to low temperature, has long 
been known to have an anomalous Mn BECfi^ sugges- 
tive of a ferroelectric instability. And second, recent 
calculation&iiii^ for CaMnOs showed that ferroelectricity 
can be artificially stabilized by increasing the lattice con- 
stants (equivalent to applying negative hydrostatic pres- 
sure) or with tensile strain, so that the short-range repul- 
sions are substantially reduced. Since the Ba^+ ion has 
a larger radius than Ca^+, perovskite-structure BaMnOa 
is analogous to CaMnOa under negative pressure. We 
note, however, that because of the large size of Ba^"*", 
the experimental structure of BaMnOa is not perovskite, 
but rather a hexagonal phase with face-shared oxygen 
octahedra. We discuss options for stabilizing the per- 
ovskite phase at the end of the manuscript. 

Our first-principles density functional calculations are 
performed within the local spin-density approximation 
(LSDA) as implemented in the Vienna ah initio simula- 
tion package (vASP)i^'2i The projector augmented wave 
method— is used with the following valence electron con- 
figurations: bs^bp^Qs^ for Ba, M^As^ for Mn and 2s'^2p'^ 
for oxygen. The Brillouin zone integrations are per- 
formed with the tetrahedron method^^ over a 9x9x9 
Monkhorst-Pack fc-point mesh^l centered at F and a 
450 eV plane-wave cutoff. For structural relaxations a 
Gaussian broadening technique of 0.05 eV is used and 
the ions are relaxed until the Hellmann-Feynman forces 
are less than 1 meV A^^. The electronic contribution 
to the polarization is calculated following the standard 
Berry phase formalism 

We begin by determining the equilibrium volume for 
the hypothetical cubic (space group Pm3m) BaMnOa 
structure with C?-type antiferromagnetic order of the 
Mn'*+ cations. This G-type antiferromagnetism is the 
likely ground state since both CaMnOa and SrMnOa 
have been shown to exhibit the same type of order 
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FIG. 1: Phonon dispersions of the G-type AFM cubic 
BaMnOa along the high-symmetry hnes. Energy lowering 
eigenvector modes for the unstable phonons at F, X and M, 
with the Mn ion at the center of the unit cell surrounded by 
an oxygen octahedron. 



in the perovskite phasei^l;^ We find that the primi- 
tive equilibrium unit cell volume fJo = 58.7 A"^ (cu- 
bic lattice parameter of — 3.89 A) is larger than 
previous first-principles calculations for CaMnOs using 
various exchange-correlation functional, consistent with 
the larger A-site cation size J^i^i^S The centrosymmetric 
structure is insulating with a LSDA band gap of 0.20 eV. 
(With a moderate effective Hubbard U of 4.5 eV, we find 
a gap of 0.26 eV) . The valence band is predominantly ma- 
jority spin Mn tig and O 2p character with strong p-d 
hybridization [Figure O (upper panels)]. The conduction 
band is formed by the empty Mn Cg orbitals and the mi- 
nority spin tig states, which is consistent with Mn'*+ in 
an octahedral crystal field. A local magnetic moment of 
2.3 /is is found at each Mn site, which is slightly reduced 
from the formal 3 /is due to strong hybridization with 
the oxygen 2p orbitals. 

We next calculate the lattice instabilities for the opti- 
mized cubic structure using the frozen-phonon method^i 
In this method the phonons are found by calculating to- 
tal energies with respect to atomic displacements from 
the reference structure at high symmetry positions, g, in 
the Brillouin zone. We use a 3x3x3 (270 atom) super- 
cell, which allows us to directly access 8 high symme- 
try points in the Brillouin zone by freezing in different 
atomic displacement patterns. For each q value we con- 
struct the dynamical matrix from the Hellman-Feynman 
forces induced on the ions after making small positive and 
negative displacements (to remove any quadratic effects) 
about the high symmetry positions. Diagonalization of 
the dynamical matrix yields the atomic displacement pat- 
terns (eigenvectors) and phonon mode frequencies (eigen- 
values) at that q point. The complete phonon dispersions 
are then determined by interpolating these solutions of 



the dynamical matrix at the special q-points to the whole 
Brillouin zone using a Fourier interpolation schemci^i^ 
The splitting between the longitudinal optic (LO) and 
transverse optic (TO) modes at fc = is not included in 
our calculations. 

In Figure [TJ we show the phonon dispersion curves for 
G-type BaMnOs along the high-symmetry directions in 
the Brillouin zone of the primitive (5-atom) unit cell. 
The dominant instability is a triply degenerate F-point 
mode with an imaginary frequency of 157. 2i cm~^ and 
point symmetry Ti„. This is a polar mode consisting of 
a relative Mn-0 displacement in which the oxygen octa- 
hedra remain almost rigid. In the ground state structure 
(described below) we find that this mode dominates with 
the Mn displacement occurring along the [110] direction 
of the primitive cubic unit cell. Interestingly, the other 
high symmetry instabilities - at the X and M points - 
do not correspond to the usual antiferrodistortive rota- 
tions of the oxygen octahedra which are common in per- 
ovskites. The X-point mode has an imaginary frequency 
of 125. 5i cm~^ and Aiu symmetry. It corresponds to 
a breathing of the oxygen octahedra in the a;?;-plane in 
which two adjacent oxygens move in towards the Mn ion 
and two move out, so that two short and two long Mn-0 
bonds are created (Figure [J). In the ground state struc- 
ture it combines with the Ti„ mode with almost equal 
amplitude. The less unstable M-point [w = 19.64z, (iJ„ 
symmetry)] also corresponds to a distortion rather than 
a rotation of the oxygen octahedra, this time dominated 
by displacements of the apical oxygen atoms (Figure [1]). 
We also verified that these modes are robust to correla- 
tion effects by repeating the calculations with an effective 
Hubbard U parameter of 4.5 eV. 

Next we determined the ground state structure by 
freezing linear combinations of each of the unstable 
modes described above into the cubic Pm3m phase, 
then fully relaxing all internal degrees of freedom un- 
til the forces were less than 1 meV A~^. We began 
with a constant volume constraint, and obtained a po- 
lar structure with R3m space group. '^^ We then lifted 
the volume constraint and performed a full optimization 
of the atomic positions and lattice parameters within 
orthorhombic symmetry (20-atom unit cell). We ob- 
tained a polar ground state structure with Amm2 space 
group, and structural parameters given in Table [D The 
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-0.011 


0(2) 


4e 


1 

2 


0.249 


-0.260 



TABLE I: Calculated structural parameters for BaMnOa with 
Amm2 symmetry. Our calculated lattice parameters are are 
a — 3.84, b — 5.43, and c = 5.43 A with a volume per formula 
unit of 56.7 A^ 
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FIG. 2: Orbital resolved Mn 3d and O 2p densities of states in 
the centrosymmetric cubic (PmSm) and the polar {Amm2) 
structures. (Upper panels) spin-polarized DOS of a single Mn 
atom in the G-type AFM BaMnOa structure. (Lower panels) 
DOS of the O atom on top of the Mn atom. In the Pm3m 
structure only the pz orbital is shown as the oxygen 2p orbitals 
are degenerate. 



ground state structure is 64 meV lower in energy per 
formula unit than the cubic Pm3m structure. We find 
that our calculated ground state structure can be writ- 
ten as a linear combination of the unstable modes of the 
cubic reference structure with the following coefficients: 
0.694Ti„ + 0.604A2„ + 0.393£'„. The resulting structure 
has three unique Mn-0 bond lengths of 1.88, 1.92 and 
1.96 A with a mean bond length similar to the cubic 
case. Furthermore, the mean Mn-O-Mn bond angle is 
decreased to 178.2° from the ideal 180° cubic structure. 

Next, we examine the changes in the electronic struc- 
ture following displacement in order to explain the stabi- 
lization of the ferroelectric off-centering of the magnetic 
Mn cation. In Figure[2]we compare the densities of states 
(DOS) for the non-polar cubic and polar Amm2 struc- 
tures; in particular, we focus on the orbital resolved DOS 
for the Mn 2>d and O 2p bands. The electronic struc- 
ture of the metastable cubic phase was described earlier; 
here we reiterate that both the t2g and eg states are hy- 
bridized with the O 2p orbitals indicating that some cova- 
lent bonding is already present in the cubic phase. In the 
polar Amm2 ground state an increase in hybridization 
between the Mn Cg and O 2p levels occurs. This lowers 
the energy of the eg states by ^0.30 eV, and shifts the O 
Pz states from the top to the bottom of the valence band; 
both factors result in an increase of the valence band- 
width. In terms of our earlier discussion of the SOJT 
effect, the matrix elements (0| TiS^^ \n) are large because 
of the enhanced Mn eg and O 2p hybridization, and the 
ferroelectric distortion is favored. 



We next calculate the total polarization for the 
Amm2 structure as the sum of the ionic and elec- 
tronic contributions (including both spin channels) us- 
ing the Berry's phase method. Here we find a value of 
12.8 /iC cm~^, which is substantial compared to many 
manganite multiferroics ,'^^'^^'^^ but consistent in magni- 
tude with conventional cP ferroelectric perovskite oxides. 
We next calculate the BECs for the cubic phase along 
the [111] direction, since large anomalies from the formal 
charge values are often good indicators for the under- 
lying ferroelectric instability. In Tabic |TT] wc tabulate 
the spin-resolved BECs for each atom i decomposed such 
that Z* = Z]™ + Zf'^ + Zf'^, where Z]™ is the pseudo- 
core charge. We find an anomalously large BEG for Mn, 
consistent with our finding of a ferroelectric instability 
dominated by Mn displacement. In fact we see that the 
majority spin electronic contribution is positive, corre- 
sponding to a net flow of electrons towards the Mn as it 
displaces towards the oxygen; this is consistent with the 
enhanced Mn-O hybridization that we observed in our 
calculated DOSs. [An anomalously large value is also 
found for Mn (8.45), when correlation effects are added, 
albeit reduced due to the enhanced band gap and nar- 
rowing of the eg bandwidth.] Interestingly, the minority 
spin electronic contribution on the Mn ion is close to the 
formal charge value; this reflects the low availability of 
Mn minority states to accept electrons at the bottom of 
the conduction band in the cubic structure. In the po- 
lar Amm2 structure, the calculated Mn Z* in the [111] 
Cartesian direction is 8.34 (7.45 within the LSDA-|-f7 
method); this reduction from the cubic case is consis- 
tent with behavior in conventional perovskites and also 
reflects in this case the opening of the band gap between 
the cubic and ground state structures. 

We conclude by exploring the feasibility of accessing 
perovskite-structure BaMnOa experimentally. As we dis- 
cussed above, non-ferroelectric CaMnOa exists in the per- 
ovskite structure and is nearly cubio^^ with small ro- 
tations of ideal MnOg octahedra. Increasing the size 
of the ^-site cation, which we have shown promotes 
the tendency to ferroelectric instability, also destabilizes 
the perovskite structure. SrMnOa is found experimcn- 
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-3.97 


6 


-4.99 


-4.98 



TABLE 11: Spin-resolved Born effective charges for cubic 
BaMnOa calculated within the LSDA along the [111] direc- 
tion. The spin components for Ba and O are equivalent within 
the numerical accuracy. 
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FIG. 3: Total energies for the the polar Amm2 and experi- 
mental 4H BaMnOa structures as a function of the pseudo- 
cubic lattice parameter. The lines are quadratic fits to the 
calculated energies. 

tally in both cubic perovskite and the non-perovskite 
2H hexagonal polymorphs?^^ while BaMnOs is known 
only in the denser hexagonal 4H structure with face- 
shared octahedra.^° Therefore, we next calculate the rela- 
tive stability of our predicted Amm2 perovskite-structure 
BaMnOs and the experimental 4H structure. Our re- 
sults are shown in Figure [31 as a function of pseudo- 
cubic lattice parameter. We find that the experimental 
4H-BaMn03 structure is more stable than the Amm2 
perovskite structure by ~0.50 eV per formula unit. Only 
near 4.7% expansion of the 4H lattice parameter is the 
perovskite structure energetically favored over the 4H 
structure. These large energy differences suggest that 
experimental synthesis of the perovskite phase will be 
challenging, although might be feasible through epitaxial 
stabilization in ultra-thin films; this approach has indeed 
been successful in stabilizing orthorhombic, over hexago- 
nal, rare-earth manganites 4 - -i -' ? ' Another possibility is to 
search for a critical alloying range with Ca and/or Sr at 
which the perovskite structure remains stable over the 
2H or 4H structure, but within the perovskite structure 



where the ferroelectric instability has already started to 
dominate. 

Finally, we note that, if a ferroelectric Ba-based per- 
ovskite manganite could be achieved experimentally, sub- 
sequent rare-earth doping on the A-site could lead to 
an intriguing analogy to the colossal magnetoresistive 
(CMR) materials. In conventional CMR materials, a 
phase transition between a paramagnetic insulator and a 
ferromagnetic metal is induced by applied magnetic field 
and/or temperature. In rare-earth doped BaMnOa, the 
insulating phase might be expected to be polar, leading 
to a ferroelectric insulator - ferromagnetic metal transi- 
tion. We suggest Ba-rich Lai-^jBa^^MnOa as a material 
for further experimental investigation. 

To summarize, we have used hypothetical perovskite- 
structure BaMnOa to demonstrate that, in contrast to 
conventional wisdom, non- magnetic cations can un- 
dergo SOJT off-centering distortions resulting in polar 
ground states. We have explained the behavior in terms 
of the delicate balance between competing energies in the 
second-order Jahn- Teller effect, and shown that reducing 
Coulombic repulsions via increasing ionic separations is 
a general route to promoting ferroelectricity in magnetic 
compounds. Finally, we have explored the feasibility of 
realizing such novel single-phase multiferroics experimen- 
tally. 
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